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D ■ Quantum corrections to the transverse spin-fluctuation propagator are obtained by 

, including self-energy and vertex corrections to first order within a spin-rotationally- 

symmetric inverse-degeneracy (1/J\F) expansion scheme which preserves the Gold- 
stone mode order by order. A correlation-induced exchange-energy correction is 
shown to yield the dominant quantum reduction in the spin stiffness, providing a 

o. 

O , quantitative understanding of the ferromagnetic-state stability in terms of simple 

lattice-dependent features of energy-band dispersion. The quantum reduction factor 
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\Q , between interaction U and bandwidth W. 
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I. INTRODUCTION 

The continuous spin-rotation symmetry of a magnetic system is manifested, in the 
spontaneously broken-symmetry state, in gapless spin-wave excitations in accordance with 
the Goldstone theorem,- the zero-energy infinite-wavelength mode simply corresponding 
to a uniform rotation of all spins. The gapless spin-wave spectrum has particularly im- 
portant consequences for low- dimensional (D = 1, 2) magnetic systems, where transverse 
spin fluctuations diverge at any finite temperature, resulting in absence of long-range 
magnetic order in accordance with the Mermin- Wagner theorem,- and exponentially large 
spin-correlation length in two dimensions. 

A quantitative determination of the spin-wave spectrum also allows for various order- 
ing, dimension, and lattice-specific investigations — finite-temperature spin dynamics and 
reduction of magnetic order with temperature due to thermal excitation of spin waves,-^ 
estimation of the magnetic transition temperature T c from the broken-symmetry side,- 1 ^ 
competing interactions and magnetic stability, 8 ' 9 and dispersion of magnetic and elec- 
tronic excitations in solids as inferred from inelastic neutron-scattering and angle-resolved 
photoemission studies,— dLlZMM References given above illustrate recent spin-wave ap- 
plications to the nearly square- and triangular-lattice antiferromagnets such as cuprates 
(La 2 Cu04), multiferroics (HoMnOa), and organic systems k — (BEDT — TTF) 2 X, as 
well as ferromagnets (Fe,Ni) and magnetic multilayers (Fe/Cr) exhibiting giant magneto- 
resistance. 

For band ferromagnets such as Fe and Ni, there have been extensive inelastic neutron- 
scattering studies in relation with calculations for transverse spin fluctuations in the ran- 
dom phase approximation (RPA), which is the lowest-order treatment in which spin ro- 
tation symmetry and Goldstone mode are preserved. Various spin-wave features such as 
isotropy, stiffness constant, damping and disappearance at higher energy due to interaction 
with the continuum of Stoner excitations, persistence for T > T c , temperature dependence 
of dispersion etc. have been discussed,— and also quantitatively compared with 
RPA calculations using realistic band structure.— 1 ^ However, despite the extensive study 
of magnetic excitations in metals and alloys over the years,— a spin-rotationally-symmetric 
extension of RPA including self-energy and vertex corrections has not been carried out 
quantitatively for fcc-type lattices with respect to ferromagnetic-state stability, spin-wave 
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and Stoner excitations, damping, and transition temperature. 

Some of the related developments beyond RPA are summarized below. The effect of 
large spin fluctuations in a nearly ferromagnetic Fermi liquid has been studied in the 
context of spin waves in He 3 within an extension of the paramagnon model beyond RPA, 
where ambiguities of the paramagnon model were shown to be resolved.— In the context 
of self-energy corrections in a band ferromagnet, the importance of vertex corrections in 
restoring the spin-rotation symmetry and Goldstone mode has been recognized at a formal 
level, and a Ward identity connecting vertex corrections to self-energy corrections has been 
derived.— While spin-wave excitations for arbitrary wave vector were not quantitatively 
discussed, the spin-wave stiffness constant was shown to be reduced from its RPA value, 
and also compared with earlier studies 2 ^ 1 ^ in the context of stability of the ferromagnetic 
state.— A variational approach has been used to improve the RPA result for the energy of 
long wave length spin- wave modes. 28 A spin- wave damping term proportional to q 6 due to 
scattering off particle-hole excitations has been obtained for a parabolic band.~ Recently 
self-energy corrections have been incorporated in a modified RPA approach, although the 
q, cu-dependence of vertex corrections was not included.- 

In this paper we provide a concrete extension beyond RPA for the transverse spin- 
fluctuation propagator. We make use of the inverse-degeneracy (1/W) expansion within 
the generalized A^-orbital Hubbard model,— which provides a systematic diagrammatic 
scheme for incorporating quantum corrections while preserving spin-rotation symmetry 
and hence the Goldstone mode order-by-order. This spin-rotationally-symmetric scheme 
has been applied earlier to examine quantum corrections in the antiferromagnetic state 
of the Hubbard model.— The diagrams include self-energy and vertex corrections, and 
physically incorporate effects such as quasiparticle damping, spectral-weight transfer and 
coupling of spin and charge fluctuations. We consider the special case of a saturated band 
ferromagnet, in which the absence of minority-spin particle- hole fluctuations results in 
relative simplification. 

Owing to its intrinsically strong-coupling nature, band ferromagnetism has been recog- 
nized as a fairly challenging problem, particularly with respect to the estimation of Curie 
temperature for the Hubbard model, although considerable progress has been achieved 
in the recent past. 31 Competition between band and interaction energies, separation of 
moment-melting and moment-disordering temperature scales due to strong correlation, 
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and presence of charge fluctuations even in the broken-symmetry state due to partially 
filled band(s) are some of the non-trivial elements involved. Ferromagnetism in the Hub- 
bard model on fee and bee lattices has been recently investigated using several different 
approaches, such as the dynamical mean field theory (DMFT),— by incorporating spin and 
charge fluctuations in the correlated paramagnet using the fluctuation-exchange (FLEX) 
and the two-particle self consistent (TPSC) approximations, 33 by systematically improv- 
ing the self energy— and a modified RPA scheme.- Ferromagnetism in a diluted Hubbard 
model has also been investigated recently,— 1 ^ which is of interest in the context of carrier- 
mediated ferromagnetism in diluted magnetic semiconductors such as Gai_ x Mn x As. 

Incorporating only the local (Ising) spin excitations, the DMFT approach ignores long- 
wavelength spin fluctuations and the k-dependence of self energy. FLEX incorporates 
self-energy corrections, but ignores vertex corrections of the same order, thereby breaking 
the spin-rotation symmetry. Both DMFT and FLEX are hence not in accordance with 
the Mermin- Wagner theorem. While self-energy corrections in the broken-symmetry state 
were incorporated in the modified RPA approach,- the momentum-energy dependence of 
vertex corrections was not considered. While FLEX, DMFT, and RPA results for the 
behaviour of Curie temperature with band filling are found to be qualitatively similar, 
appreciable quantitative differences^ 13 ^ 3 - clearly highlight the need for a spin-rot at ionally- 
symmetric extension. 

We consider the generalized A/"-orbital Hubbard model 3 ^ 

H = ~t (4m a j^ + H - c -) + jj^{ U ^ a \\a a iU a \ii3 a iW + U 2 a\ }a a iW a\ ip a^ a ) , (1) 

(ij),a,a i,o,f3 

where a, (3 refer to the degenerate orbital indices and the factor 1/jV is included to render 
the energy density finite in the M — ► oo limit. In the isotropic limit U\ = U 2 = U, the 
two interaction terms (density- density and exchange-type with respect to orbital indices) 
are together equal to U(— Sj.Sj + nf) in terms of the total spin S« = J2 a i J la(, <T /'2)ipia 
and charge rii = ^2 a ip} a (l/2)il>ioi operators, and the Hamiltonian is therefore explicitly 
spin-rotationally symmetry. 
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II. TRANSVERSE SPIN FLUCTUATIONS 

The transverse spin-fluctuation propagator in the broken-symmetry state, which de- 
scribes both collective spin-wave and particle-hole Stoner excitations, is given by 

X -+(q,c) = i fdt e^J^E^'^^lT^W-SjCOl^a) (2) 

in terms of the fermion spin-lowering and -raising operators S T = ^ t (cr =F /2)^. The 
spin-fluctuation propagator can be expressed as 

1 - U(p{q,uj) 

in terms of the exact irreducible propagator </>(q, uj), which incorporates all self-energy 
and vertex corrections. The inverse- degeneracy expansion 30 

2 # (2) + ... (4) 



systematizes the diagrams in powers of the expansion parameter 1/JV which, in analogy 
with 1/S for quantum spin systems, plays the role of h. As only the "classical" term 
survives in the Af — > oo limit, the RPA ladder series x°(q, cj)/1 — f/%°(q, uj) (with 
interaction C/2) amounts to a classical- level description of non-interacting spin-fluctuation 
modes. The bare antiparallel-spin particle-hole propagator 

0(°)(q,c) = X °(q,c) = V— i , (5) 

k 4- q -4 +^-^ 

where = ek — crA are the Hartree-Fock ferromagnetic band energies, 2A = mU is the 
exchange band splitting, and the superscript +(— ) refer to particle (hole) states above 
(below) the Fermi energy ep. For the saturated ferromagnet, the magnetization m is equal 
to the particle density n. 

As collective spin- wave excitations are represented by poles in (3), spin-rotation sym- 
metry requires that = 1/U for q, uo = 0, corresponding to the Goldstone mode. Since the 
zeroth-order term <p°) already yields exactly 1/U for q, uj = 0, the sum of the remaining 
terms must exactly vanish in order to preserve the Goldstone mode. For this cancel- 
lation to hold for arbitrary J\f, each higher-order term cf)^ in the expansion (4) must 
individually vanish, implying that spin-rotation symmetry is preserved order-by-order, as 
expected from the spin-rotationally-invariant form (U/J\f)Si.Si of the interaction term 
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FIG. 1: The first-order quantum corrections to the irreducible particle- hole propagator 0(q, w). 

in the generalized Hubbard model.— We evaluate the order 1/jV diagrams in 0^ and 
explicitly show the exact cancellation for q, uo = 0. 

We consider the relatively simpler case of a saturated ferromagnet in which the 
minority-spin (|) band is pushed above the Fermi energy due to Coulomb repulsion, 
resulting in the absence of any minority-spin particle-hole processes. In this case the 
effective antiparallel-spin interaction at order 1/Af reduces to the bare Hubbard interac- 
tion U and the effective parallel-spin interaction reduces to a single term involving the 
majority-spin (f) particle-hole bubble. Generally these effective interactions involve a 
series of bubble diagrams, with even and odd number of bubbles, respectively. 

The order 1/M diagrams for the irreducible particle-hole propagator </>(q, u) are 
shown in Figure 1. The hatched part in diagram (a) represents the RPA ladder sum 
^ 2 Xrpa(~Q> ~ ^) with interaction U 2 , where 

x°(Q,n) 



i-u x °(Q,n) 



(6) 
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has purely advanced (retarded) character with respect to Q (— Q), and includes both spin- 
wave and Stoner excitations. Diagram (a) involves a self-energy correction to the j-spin 
particle which transfers spectral weight down from the j-spin band to the j-spin band 
(above ep), and yields a positive correction to <fi. Diagrams (b) and (c) represent vertex 
corrections, where the hatched part represents U/l — f/x°(Q,f2), the RPA ladder series 
starting with a single interaction line U. In diagram (b) the single opposite-spin particle- 
particle interaction (Ui) reduces the j-spin particle — j-spin hole correlation, yielding a 
negative correction to 0, whereas two such {Ui) interactions in (c) and (d) yield positive 
corrections. 

A coupling between spin and charge fluctuations is indicated by the j-spin particle-hole 
bubble, present explicitly in diagrams (c) and (d) and implicitly in (a) and (b). It is the 
availability of unoccupied j-spin states for partial band filling n < 1 which allows for the 
different processes (a)-(d), either involving spin flip due to spin- wave coupling (a) or the 
fermion-fermion scattering due to on-site Coulomb interaction (b,c,d). Indeed all these 
quantum corrections identically vanish for a completely filled j-spin band. This coupling 
between spin and charge fluctuations is a unique feature of the quantum corrections and 
is absent at the classical level. 

Integrating out the fermion frequency-momentum modes, the order 1/M quantum 
corrections to the irreducible particle-hole propagator 0(q, uj) are obtained as: 

x°(Q,n) 



(a) (q,^) 

Q 




tfx°(Q,n) 
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In the infinite-wavelength limit (g — > 0), the total order l/J\f contribution ^^(q, uj) 
exactly vanishes, as required from spin-rotation symmetry. With e k + — e k ~ = 2A, we 
obtain 

{1) (q = 0,u) = (a) + + {c) + 



= E rxi/ i \w 



(11) 



_l-[/ x o 1 _ Ux o 1 _ Ux o 

which yields identically vanishing contribution for each spin-fluctuation mode Q. We 
note that this mode-by-mode cancellation is quite independent of the spectral- weight 
distribution of the spin-fluctuation spectrum between collective spin-wave excitations and 
particle-hole Stoner excitations. Furthermore, the cancellation holds for all uj, indicating 
no spin-wave amplitude renormalization, as expected for the saturated ferromagnet in 
which there are no quantum corrections to magnetization. 

The coupling between spin and charge fluctuations is highlighted by the structure 
of (11), where the common f-spin particle-hole bubble term within the k' sum repre- 
sents charge fluctuations. This coupling provides a spin-wave damping mechanism for a 
band ferromagnet, fundamentally different from the conventional damping mechanism in 
Heisenberg insulating magnets involving decay into three spin waves. The spin-wave de- 
cay and damping process in a band ferromagnet involves the imaginary part of the j-spin 
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particle-hole bubble 

EE^ + Q- £ t+ w - n ) ( 12 ) 

Q k' 

which corresponds to energy conservation in the spin-wave (energy — uj = uj-q) decay 
into an intermediate-state spin wave (energy —Q = ^-q) plus a particle-hole excitation 
in the metal. Exactly vanishing for uj = 0, the imaginary term increases with |u;|. The 
typically strongly peaked spin-wave density of states near the top end of the spectrum 
yields significant spin-wave damping only for zone boundary modes. Further investiga- 
tions, including effects of disorder and diffusion pole, are clearly of interest in view of the 
observed temperature-independent linewidth in neutron-scattering studies.— 1 ^ 



III. QUANTUM CORRECTIONS TO SPIN STIFFNESS 

We next consider the quantum corrections for small q and show that the above exact 
cancellation for q = actually extends to the next order as well, yielding spin-stiffness 
quantum corrections only from the higher-order surviving terms. For analytical simplicity, 
we neglect the contribution of higher-energy Stoner excitations in the Q integral in (7) - 
(9), and only consider contributions from the spin- wave pole 

X°(Q,n) rn Q 

i-u x °(Q,n) Q + Q Q -V 1 ' 

where rriQ is the spin-wave amplitude for momentum Q. Writing the antiparallel-spin 
particle-hole energy denominators as 

4- q - 4~ = 2A[1 + (e k _ q - e k )/2A] (14) 
and expanding in powers of the small band-energy difference 

5 = -(e k _ q - e k ) = q.Ve k - ^(q.V) 2 e k (15) 

for small q, we find that besides the zeroth-order cancellation for q = 0, the first-order 
terms in 5 also exactly cancel, implying no quantum correction to the classical term 
(V 2 e k ) in the spin-wave stiffness constant, in accordance with the exact structure.— 1 ^ 
The surviving second-order terms in 5 can be written, up to order q 2 , as 



U* (q.Ve k 2 
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where we have set q, uo = in the energy denominators as all three terms are already 
explicitly second order in q. Incorporating the energy cost of spin twisting, the three 
terms in (16) represent exchange-type processes involving f-spin particle-hole (charge) 
excitations accompanied with zero, one, and two j-spin particle-particle scatterings, re- 
spectively. Cross terms such as q x q y etc. in (16) identically vanish from symmetry, leaving 
an isotropic momentum dependence on q 2 = g 2 + q 2 + q\. Equation (16) is a new result 
and incorporates all first-order (1/JV) quantum corrections to the spin-wave stiffness con- 
stant due to collective spin-wave excitations. Including the Stoner contribution in (7-10) 
again yields an exact cancellation of the V 2 €k-type terms, leaving only an additional 
(Vtk) 2 -type term qualitatively similar to the first term in (16). 36 

Now, as Vek is odd in momentum k, the second and third terms in (16) involve a 
partial cancellation resulting from the momentum summations. Therefore, keeping the 
contribution of the dominant first term only, with an order-of-magnitude estimate for 
the particle-hole energy denominator as W, the fermion bandwidth, the renormalized 
spin-wave energy in D dimensions is obtained as 



2 V k/ 2A 2A W y 1 



Q 2 , (17) 



where 2A = nil, the angular bracket ( ) represents momentum summation normalized 
over the number of occupied states, the explicit hole density factor (1 — n) highlights the 
particle-hole process involved which vanishes for a filled band, a is a band-dependent factor 
of order 1, and the prefactor 1/D follows from hypercubic symmetry. A straightforward 
calculation of the two competing terms in spin stiffness (V 2 e k ) and ((Ve k ) 2 )/2A can 
provide a quantitative estimate of the stability of the ferromagnetic state for different 
lattices. 

The origin and physical interpretation of the three terms in (16) are given below. 
The first two terms represent classical (RPA) spin-stiffness contributions arising from 
the two fermion band curvature terms (15) substituted in (5). With respect to stiffness 
against spin twisting, the first (positive) term represents delocalization-energy loss which 
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vanishes for a filled band, whereas the second (negative) term of order t 2 /U represents 
exchange-energy gain. The third term represents an additional exchange process involving 
minority-spin intermediate states which are transferred to lower energies, corresponding 
to the finite probability 1 — n of a site being unoccupied by a majority-spin electrons. The 
third term thus represents a correlation-induced quantum correction to spin stiffness. 

If the normalized averages () in (17) are essentially n-independent (as for the fee lat- 
tice with t' = 0.25 in the low-density limit), the band-filling dependence (1 — n)/n of the 
quantum term destabilizes the ferromagnetic state for small n, whereas the competition 
between the two classical terms destabilizes the ferromagnetic state as n approaches 1, 
yielding an optimization of the spin stiffness at some intermediate n. Also, the quantum 
correction factor U/W again highlights the characteristic competition between interaction 
U and bandwidth W, although favouring a stability condition quite opposite to the Stoner 
criterion. In the low-density limit (n«l), the second and third terms in (17) can be com- 
bined as ((Vek) 2 ) /nU e s, where the effective interaction U e s = U/(l + U/W) approaches 
the bandwidth W in the strong coupling limit, in agreement with the low-density result 
of Kanamori.— 

IV. CONCLUSIONS 

We have investigated correlation effects on the spin dynamics in a band ferromagnet, 
providing a physically transparent framework for quantitative understanding of ferromag- 
netic stability in terms of simple lattice-dependent features of energy-band dispersion. 
The correlation effect arises from the minority-spin spectral-weight transfer to lower en- 
ergies corresponding to finite site- vacancy probability (1 — n) of majority-spin electrons, 
the availability of these low-lying intermediate states resulting in an additional exchange- 
energy gain of order t 2 /U, yielding a lattice-dependent destabilization of the ferromagnetic 
state. 

The inverse-degeneracy (1/W) expansion provides a fully spin-rotationally-symmetric 
scheme for incorporating self-energy and vertex corrections in the transverse spin- 
fluctuation propagator. First-order contributions to the irreducible particle-hole prop- 
agator were obtained, with full momentum-energy dependence in the vertex corrections, 
and shown to have appropriate cancellations for q = 0, small q, as well as finite us. 
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Lattice-specific evaluations thus allow for quantitative study of magnetic excitations for 
arbitrary wave vector, which should be especially suitable for low-dimensional systems as 
the scheme is in accordance with the Mermin- Wagner theorem. 

The factor U/W in the spin-stiffness quantum reduction highlights the subtlety in the 
characteristic competition in a band ferromagnet between interaction U and bandwidth 
W. While Stoner criterion favours a high density of states at Fermi energy, the quan- 
tum correction favours a large bandwidth for the stability of the ferromagnetic state. 
Indeed, with a large density of states at one end and a broad band tail, the fcc-lattice 
band does provide optimum conditions for a stable ferromagnetic state within the single- 
band Hubbard model. The quantum corrections also involve a coupling between spin 
and charge fluctuations, resulting in an intrinsic spin-wave damping mechanism which 
fundamentally distinguishes between band and insulating ferromagnets. Lattice-specific 
calculations including the contribution of both spin-wave and Stoner excitations in the 
intermediate-state spin-fluctuation spectrum are currently in progress, and quantitative 
results for spin stiffness, spin-wave energy, Curie temperature, electronic spectral function 
etc. will be presented elsewhere.— 
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